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Abstract 

Background: The yeast Xanthophyllomyces dendrorhous is used for the microbiological production of the 
antioxidant carotenoid astaxanthin. In this study, we established an optimal protocol for protein extraction and 
performed the first proteomic analysis of the strain ATCC 24230. Protein profiles before and during the induction 
of carotenogenesis were determined by two-dimensional polyacrylamide gel electrophoresis and proteins were 
identified by mass spectrometry. 

Results: Among the approximately 600 observed protein spots, 131 non-redundant proteins were identified. 
Proteomic analyses allowed us to identify 50 differentially expressed proteins that fall into several classes with 
distinct expression patterns. These analyses demonstrated that enzymes related to acetyl-CoA synthesis were more 
abundant prior to carotenogenesis. Later, redox- and stress-related proteins were up-regulated during the 
induction of carotenogenesis. For the carotenoid biosynthetic enzymes mevalonate kinase and phytoene/squalene 
synthase, we observed higher abundance during induction and/or accumulation of carotenoids. In addition, 
classical antioxidant enzymes, such as catalase, glutathione peroxidase and the cytosolic superoxide dismutases, 
were not identified. 

Conclusions: Our results provide an overview of potentially important carotenogenesis-related proteins, among 
which are proteins involved in carbohydrate and lipid biosynthetic pathways as well as several redox- and stress- 
related proteins. In addition, these results might indicate that X. dendrorhous accumulates astaxanthin under 
aerobic conditions to scavenge the reactive oxygen species (ROS) generated during metabolism. 



Background 

The basidiomycete Xanthophyllomyces dendrorhous (for- 
merly known as Phaffia rhodozyma) is an excellent 
astaxanthin-producing yeast and has been regarded as 
one of the most promising microorganisms for the com- 
mercial production of this carotenoid [1,2]. Astaxanthin 
is a pigment that produces the characteristic coloration 
of some birds, crustaceans and salmon. It has been used 
as a feed and food pigment in the aquaculture industry 
and has been evaluated as a pharmaceutical component 
because it may possess antioxidant activity [3,4]. Due to 
its biotechnological significance, investigations have 
been performed to improve astaxanthin production by 
optimizing fermentation methodologies [5,6] selecting 
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for over-producing strains [7,8], using chemical stimu- 
lants [9,10], and employing genetic and metabolic engi- 
neering [11-13]. 

In X. dendrorhous, astaxanthin is produced via the 
mevalonate pathway, in which acetyl-CoA is a precursor 
to the formation of isopentenyl pyrophosphate (IPP), the 
general precursor of all isoprenoids. Chain elongation by 
successive head-to-tail condensation is then catalyzed 
by prenyltransferases with different chain length specifi- 
cities [14] to generate carotenoid precursors. Thus, 
farnesyl pyrophosphate (FPP) (CI 5) and geranyl pyro- 
phosphate (GGPP) (C20) are the immediate precursors 
of C30 and C40 carotenoids. GGPP formation is cata- 
lyzed by a GGPP synthase. The condensation of two 
molecules of GGPP is catalyzed by the bifunctional 
enzyme phytoene synthase to produce phytoene (C40). 
Lycopene is generated by phytoene desaturase, which 
introduces four double bonds into phytoene. A 
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bifunctional enzyme with a lycopene cyclase activity 
then transforms lycopene into P -carotene by two cycli- 
zation reactions. Finally, P -carotene is oxidized by astax- 
anthin synthetase to yield astaxanthin [15]. 

Because little information on the genomics and regu- 
lation of carotenogenesis in X. dendrorhous is available, 
studies of astaxanthin production from a genetic per- 
spective have been hampered; however, an alternative 
approach to address these biological questions is proteo- 
mics Two-dimensional (2D) techniques are the most 
generally applicable methods for obtaining a global pic- 
ture of protein expression levels, and mass spectrometry 
(MS) has become the technology of choice for protein 
identification [16,17]. In previous studies, it has been 
demonstrated that 2D electrophoresis coupled with pep- 
tide mass fingerprinting (PMF) is a viable approach for 
the identification of homologous proteins across species 
boundaries [18-21]. Therefore, biosynthetic pathways 
and metabolic events in X. dendrorhous may be deduced 
from the functions of previously identified proteins. 

In the present study, we used 2D protein electrophor- 
esis coupled with matrix-assisted-laser-desorption/ioniza- 
tion time-of-flight mass spectrometry (MALDI-TOF MS) 
to analyze soluble protein extracts from X. dendrorhous 
cells grown on glucose minimal medium (MM-glucose). 
To the best of our knowledge, this is the first proteomic 
study on this yeast; thus, prior to protein characteriza- 
tion, we designed an optimized protocol for protein 
extraction. Because some specific or late reactions in car- 
otenogenesis involve membrane-bound enzymes, we 
designed a protocol for the enrichment of membrane- 
bound proteins. These extracts were separated in two 
dimensions to obtain a protein map. In our analysis, the 
most abundant proteins were involved in primary meta- 
bolic pathways, and carbohydrate and lipid metabolic 
proteins showed the highest intensity spots. Interestingly, 
along with some carotenogenesis proteins, redox- and 
stress-associated proteins were up-regulated. This pro- 
teomic study is an important starting point and may be a 
useful reference for further studies of metabolic path- 
ways, especially astaxanthin synthesis in X. dendrorhous. 

Results and discussion 

Isolation of soluble proteins and 2D electrophoresis 

The aim of the present study was to characterize the 
proteome of soluble protein extracts of the yeast X. den- 
drorhous grown in MM-glucose media. Therefore, sam- 
ple preparation prior to the IEF run was important for 
optimal electrophoretic separation of proteins. 

In the Methods, we describe the comprehensive proto- 
col used to obtain the soluble protein extracts. Briefly, 
to improve cell disruption and minimize proteolysis, lyo- 
philized yeast cells were vortexed directly with glass 
beads. Lysis buffer and protease inhibitors were then 



added to reduce proteolytic enzyme activity. The pellet 
was disrupted five times in a RiboLyzer, followed by 
phenol extraction and methanol precipitation. Finally, 
the protein spots were stained with Coomassie and 
identified by MALDI-TOF MS. 

To obtain the protein profiles of X. dendrorhous, the 
yeast was cultured in MM-glucose and harvested at the 
lag, late exponential and stationary growth phases. Four 
independent cultures showed continuous increases in 
cell density until 70 h, which was immediately prior to 
the induction of carotenoid biosynthesis (Figure 1). As 
we have previously reported, pigment accumulation in 
MM-glucose was evident during the stationary phase 
[22,23]. Carotenoid analysis by HPLC showed that astax- 
anthin was the main carotenoid (75-90% of the total 
carotenoids) produced by the yeast during growth. 

For the proteomic analyses, triplicate protein extracts 
(prepared from three independent cultures) were sub- 
jected to 2D analysis, and their protein profiles were 
obtained. The different protein profiles were subjected 
to a stringent comparative analysis using PDQuest soft- 
ware (version 7.1.1, Bio-Rad). After automated spot 
detection, spots were checked manually to eliminate 
possible artifacts such as background noise or streaks. 
Student's t-test (p < 0.05) was used to determine 
whether the relative changes in protein abundance were 
statistically significant. 

A representative 2D image is shown in Figure 2. The 
protein data analyses showed a consistent protein profile 




, 1 1 1 1 1 1 1 1 ] 1 1 1 1 1 1 — 0 

0 24 48 72 96 120 144 168 192 

Time (H) 

Figure 1 Growth and pigment production in X. dendrorhous 

Growth was measured by the absorbance at 560 nm (shown on a 
log scale), which is represented by the squares and solid line. The 
means ± SD of the values obtained from four independent cultures 
are shown. The vertical arrows indicate the harvest times for the 
assays (24, 70 and 96 h, which corresponded to lag, late exponential 
and stationary growth phases, respectively). The solid line represents 
the total carotenoids. The asterisk indicates the induction of 
carotenoid biosynthesis. 
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during growth (See additional file 1, Fig. SI). On aver- 
age, approximately 600 spots were detected on each 2D 
gel in a pi-range of 3-10 and a molecular mass range of 
10-100 kDa. This pattern of proteins was highly repro- 
ducible, and similar results were obtained in the tripli- 
cate cell extracts. Overall, the protein profiles did not 
change dramatically (over 90% of the spots were identi- 
cal) during growth. Of the spots detected in all gels, 450 
spots with different intensities were selected to be 
excised, digested with trypsin and analyzed by MALDI- 
TOF MS for protein identification. In total, 131 non- 
redundant proteins (we only considered one protein that 
was present in multiple spots) were identified in 171 
spots (see additional file 2, Table SI). The proteins iden- 
tified were classified according to their biological func- 
tions. Because it was impossible to determine the spot 
intensities for overlapping spots, we only quantified 161 
single-protein spots. 

Evaluation of multiple spots and differentially migrating 
proteins 

Proteins expressed from a single gene can migrate to 
multiple spots on 2D gels due to either post-translational 
modifications (such as chemical modification, proteolytic 



processing, and covalent attachment of small adducts) or 
artifactual modifications. It has been reported that several 
yeast proteins are resolved in multiple spots on 2D gels 
[24,25]. Consistent with these findings, we identified 22 
proteins that were represented by multiple spots (see 
additional file 2, Table SI), and approximately 10 pro- 
teins were present in more than three spots (Figure 2 and 
additional file 1, Fig. SI), including the stress-related pro- 
teins HSP70 (protein N°99) and ATP synthase |3 (protein 
N°82), which were previously reported to have multiple 
spots [26,27], and PP1-1 (protein N°19), a protein that 
regulates the cellular response in glucose starvation and 
stress [28]. 

In most cases, multi-spot proteins showed charge var- 
iations (horizontal spot patterns, Figure 2 and additional 
file 1, Fig. SI), which are usually due to protein phos- 
phorylation or other post translational modifications 
that alter the isoelectric point of a protein [29]. Interest- 
ingly, we found the protein, methionyl-tRNA formyl- 
transferase (protein N°69), that had a diagonal spot 
pattern, which is less frequently reported (Figure 2 and 
additional file 1, Fig. SI). This migration pattern agrees 
with the results of previous studies [24-27,30], in which 
several metabolic proteins displayed distinct migration 
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patterns. These multiple electrophoretic species could be 
generated by proteolytic events or could represent iso- 
enzymes [29]; these possibilities were not further investi- 
gated in this work. 

Approximately 25% of the proteins identified in this 
study had potential posttranslational modifications or 
belonged to multigenic protein families. Accordingly, we 
studied the intensity profiles of proteins with multiple 
spots (Figure 3), of 22 multi-spot proteins identified 8 
subgroups of proteins share similar profiles. For 
instance, a higher abundance of methionyl-tRNA for- 
myltransferase and myosin-associated protein were 
observed at the end of the exponential phase. (Figure 
3A). In contrast, acetyl-CoA carboxylase and Golgi 
transport protein showed dramatic decreases in spot 
intensity in the late exponential phase. Monooxygenase 
and kynurenine 3-monooxygenase showed increasing 
intensities during growth. Moreover, other sets of spots 
that corresponded to the same protein were notably dif- 
ferent (Figure 3B), suggesting that the isoforms are regu- 
lated in different ways or are involved in different 
physiological processes. This form of regulation has 
been previously reported for some proteins involved in 
carbohydrate metabolism [16,31]. Unfortunately, no data 
could be extracted from our MALDI-TOF analyses to 
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Figure 3 Relative intensities of multiple spots for X. 
dendrorhous proteins in MM-glucose. The growth phases are 
represented by different colors. A. Multi-spot proteins that exhibited 
essentially the same general pattern of variation. B. Multi-spot 
proteins that were regulated in different ways. The axis numbers 
correspond to the SSP spot identifications generated by PDQuest 
software. The y axis scale (x 10 3 ) corresponds to the normalized 
spot intensity. To normalize, the spot intensities were divided by the 
total density of valid spots and then multiplied by 10 6 . Finally, the 
normalized values from replicates of 24-h, 70-h and 96-h were 
averaged. Asterisks represent p < 0.01 and circles represent p < 
0.05. 



identify differences between the probable isoforms 
identified. 

Regarding the migration of proteins, for which full X. 
dendrorhous sequences were available, the experimental 
Mr and pi values corresponded closely to the theoretical 
values, except for acetyl-CoA carboxylase (N°84). For 
this protein, the experimental Mr was markedly lower 
than the theoretical Mr. This discrepancy in Mr could 
be linked to either in vivo or in vitro protein degrada- 
tion. In fact, this protein was identified with peptides 
that spanned the middle and carboxy terminal regions 
of the reported amino acid sequence. However, for the 
orthologous proteins identified, we found reasonable 
correlations between the experimental and theoretical 
migrations (see additional file 2 Table SI). Most discre- 
pancies corresponded to a lower Mr value and more 
acidic pi value for the gel-estimated value compared to 
the theoretical value. For instance, phosphatidylserine 
decarboxylase (protein N"85) was detected in the acidic 
range (pi 6.24), but this protein has a basic theoretical 
pi of 9.45. This unusual migration has been observed in 
ribosomal proteins in previous studies [30]; while this 
behavior still has no explanation, it is probably related 
to the presence of posttranslational modifications. 

Protein identification and classification into functional 
groups 

We employed the approach of cross-species protein 
identification for X. dendrorhous because this yeast is 
poorly characterized at the gene and protein levels. The 
conserved nature of many biosynthetic and metabolic 
pathways in different organisms has been the basis for 
several studies of species that lack genome sequence 
data [18,20,21]. Furthermore, based on theoretical and 
experimental evidence, it has been shown that sequence 
identities greater than 60% are sufficient to match pro- 
teins with high confidence across species boundaries 
[20,21]. It fact, it has been previously reported that con- 
served structural motifs could be identified across dis- 
tant species with total amino acid sequence identities as 
low as 29.6% [18]. In this work, nitrite reductase was 
identified with 14 mass peptides that covered 16% of the 
sequence; two of these mass peptides were located in 
the bacterioferritin-associated ferredoxin-like (BFD) 
[2Fe-2S] binding domain [32]. For mevalonate kinase, 
the four peptides identified spanned the domain desig- 
nated mevalonkin [33]. 

The proteins identified based on these analyses are 
listed in additional file 2, Table SI, along with their cor- 
responding spot numbers from the 2D gel (Figure 2). 

The proteins were classified into different groups 
according to their biological functions, which were 
determined using annotations from the KEGG database. 
The most abundant proteins found in this study were 
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involved in metabolic pathways (49%; 64 proteins) 
(Figure 4A). Others were involved in cellular transport 
(17%; 13 proteins); environmental information proces- 
sing, such as signal transduction proteins (6%; 5 pro- 
teins); genetic information processing including 
translation and transcription, replication, repair, folding 
and processing (25%; 33 proteins); and unknown pro- 
cesses (8%; 11 proteins) (Figure 4A). A similar distribu- 
tion has been observed in previous yeast proteomic 
studies (see additional file 3, Table S2). 

In the metabolism group, we identified proteins that 
belonged to different biosynthetic pathways, including 
amino acid, nucleotide, carbohydrate, energy, isoprenoid, 
redox and lipid metabolism (Figure 4B). The carbohy- 
drate-related protein group, included enzymes from the 
glycolysis, pentose phosphate (PP) and tricarboxylic acid 
(TCA) pathways (see additional file 2, Table SI). In gen- 
eral, proteins involved in carbohydrate, amino acid, 
redox and lipid metabolism showed the greatest spot 




Secondary 
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/Carotenoid 
Biosynthesis 
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Figure 4 Classification of identified proteins by cellular 
function. A. Pie chart showing the functional classifications of the 
identified proteins based on annotations from the KEGG and Swiss- 
Prot/TrEMBL protein databases. B. Proteins involved in metabolism 
(49%) were subdivided according to pathway modules in the KEGG 
database. Percentages were calculated by dividing the number of 
proteins in the group by the total number of proteins identified 
and then multiplying by 100. 



intensities when compared with all other identified 
proteins. 

Differential protein abundance during growth in MM- 
glucose 

Statistical analyses were performed using Student's t-test 
(Table 1) to select spots that showed significant changes 
in intensity relative to the intensity in the lag phase. A 
total of 66 spots (corresponding to 50 proteins) showed 
more than two-fold changes with confidence levels of 
95-99% (p < 0.05 and p < 0.01, respectively). 

Most of the differentially regulated proteins (63%) fell 
within three functional groups (metabolism, genetic 
information processing and cellular processes), while 
13% had unknown functions (Table 1). In addition, we 
observed similar patterns of intensities between proteins 
with multiple spots, such as myosin-associated protein 
and Golgi transport protein (Table 1, Figure 5). 

The fold-change data for proteins with differential 
abundances indicated that more than half of the pro- 
teins in the late exponential phase were down-regulated 
compared to their expression in the lag phase. In con- 
trast most of the proteins were up-regulated in the sta- 
tionary phase (Figure 5). Higher fold-changes were 
found for amino acid biosynthesis and transport pro- 
teins. Notably, some proteins involved in signal trans- 
duction and carotenoid biosynthesis were up-regulated 
in the late exponential and stationary phases. However, 
some redox proteins and the unknown proteins were 
down-regulated in both phases. 

In the following sections, we present an in-depth ana- 
lysis of the protein abundance patterns based on func- 
tional groups. 

Carbohydrate and lipid metabolism proteins 

In the presence of glucose, the major pathways of carbo- 
hydrate metabolism are activated to produce energy for 
the cell. Therefore, many proteins that are important for 
growing cells also play a role in stationary phase growth 
[24]. Among these proteins, the enzymes of glycolysis 
and the TCA and PP pathways were identified in the 2D 
gels. In general, this group of proteins showed high and 
similar levels of abundance during growth, which is con- 
sistent with previous reports [16,34]. As indicated in 
Figure 5 and Table 1 only two proteins (phosphogluco- 
mutase and acetyl-CoA carboxylase) were differentially 
regulated (See additional file 4, Fig. S2). It is noteworthy 
that these proteins not only have pivotal roles in central 
metabolism but are also linked to carotenogenesis. 

During the induction of carotenogenesis, phosphoglu- 
comutase (protein N°107, SSP 7519), an enzyme of the 
PP pathway, showed a three-fold increase in intensity 
(Table 1; Figure 5 and additional file 4, Fig. S2). It has 
been previously shown that astaxanthin synthesis 
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Table 1 Differentially regulated proteins of X.dendrorhous in MM glucose. 
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Table 1 Differentially regulated proteins of X.dendrorhous in MM glucose. (Continued) 
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a SSP numbers were assigned by PDQuest software analysis. 

b Identifications were obtained using the Swiss-Prot and KEGG Pathways databases and contigs of X. dendrorhous genomic DNA. 
c Data derived from PDQuest estimation. 

d Mean fold changes compared with the 24 h cultures. Bold values indicate p < 0.01, italic p < 0.02 and underlined values indicate p < 0.05. Avg., average; SD, 
standard deviation. 



requires oxygen and NADPH, which may be due to the 
reactions converting |3-carotene to astaxanthin [15]. In 
addition, the PP pathway may serve as a key source of 
NADPH for ROS removal in response to oxidative stress 
[35], and phosphoglucomutase shows changes in expres- 
sion related to NADPH generation when cells are trea- 
ted with H 2 0 2 [25]. Thus, our result suggests that high 
activity of this pathway might be required to generate 
sufficient NADPH for ROS quenching in X. 
dendrorhous. 

Acetyl-CoA carboxylase (number SSP 3516) showed a 
distinct abundance pattern during growth. This protein 
was present at high levels during the lag phase (Figure 
3A), followed by a decrease at the end of the exponen- 
tial phase and then a slight increase in the stationary 
phase. It should be noted that only one spot showed a 
significant change in intensity; the other two spots 
showed a similar trend, although these changes were 
not significant (Table 1). The decrease in abundance of 



this protein coincided with the induction of caroteno- 
genesis at the end of the exponential phase. This 
enzyme contributes to the overall control of energy 
metabolism in the cell and, catalyzes the carboxylation 
of acetyl-CoA to form malonyl-CoA, an essential sub- 
strate for fatty acid synthesis [36]. Therefore, the 
observed decrease in abundance might be related to the 
increased availability of acetyl-CoA for carotenoid 
biosynthesis. 

Although most of the carbohydrate and lipid metabo- 
lism proteins showed similar levels during growth, we 
observed that several proteins related to acetyl-CoA 
synthesis showed maximal abundance in the lag phase, 
prior to the induction of carotenogenesis (Table 1), 
including acetyl-CoA synthetase, alcohol dehydrogenase 
and ATP-citrate lyase (See additional file 4, Fig. S2) 
[37,38]. This result indicates that carbon flux to the bio- 
synthetic pathways, including carotenogenesis, is tightly 
regulated to maintain cell activity in X. dendrorhous. 
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Fold Change Fold Change 

Figure 5 Fold changes of differentially expressed proteins. Proteins with more than two-fold changes (see Table 1) were plotted according 

to their fold change in exponential phase (left graph) or stationary phase (right graph) relative to their abundance in lag phase. The proteins are 

colored by their functional classification and grouped by fold change (shown on the x-axis of each graph). SSP spot identification was 
performed using PDQuest software. 



Redox and stress response proteins 

Carotenoid accumulation is thought to be a survival 
strategy not only for the alga H. pluvialis but also for 
other microorganisms, including X. dendrorhous [39]. It 
has been observed that carotenoid biosynthesis in caro- 
tenoid-producing microorganisms is stimulated by oxi- 
dative stress [40,41]. Cellular antioxidant mechanisms 
include both non-enzymatic molecules, such as glu- 
tathione and several vitamins, and ROS scavenger 
enzymes, such as superoxide dismutase (SOD), catalase 
and glutathione peroxidase [42]. Apparently, X. dendror- 
hous lacks these enzymatic defense systems [3]; in fact, 
we identified only the mitochondrial MnSOD protein 
(see additional file 2, Table SI). This protein showed a 
higher abundance at the end of the exponential phase 
and continued to decrease during growth (Table 1 and 
additional file 4, Fig. S2). A proteomic study of H. plu- 
vialis found that this protein is constitutively highly 
expressed and is progressively down-regulated after 
stress induction (see additional file 3, Table S2). In con- 
trast, cytosolic CuSOD was found to be present in trace 
amounts and only up-regulated 48 h after stress induc- 
tion [43]. Thus, an increase in the level of CuSOD and 
modulation of the level of MnSOD were found in 
response to stress in this carotenogenic alga. Moreover, 
in a comparative analysis of C. albicans grown on glu- 
cose-supplemented media, Sod21p (cytosolic manga- 
nese-dependent) was detected only in the stationary 
phase, whereas the Sodlp isoenzyme (Cu and Zn super- 
oxide dismutase) was found only during exponential 



growth [24] (see additional file 3, Table S2). Taken 
together, these results suggest that the regulation of 
SOD is species-specific and depends on the growth 
phase. In the specific case of X. dendrorhous, we 
observed an increased level of MnSOD that coincided 
with the induction of carotenogenesis, which reinforces 
the antioxidant role of astaxanthin in the absence of 
other enzymatic antioxidant mechanisms. 

For the redox and stress response proteins, we 
observed distinct abundance patterns occurring before 
or during the induction of carotenogenesis. At the end 
of the exponential phase, three redox proteins were 
more than four-fold down-regulated (Table 1, Figure 5 
and additional file 4, Fig. S2), indicating that they were 
highly abundant in the lag phase. Interestingly, along 
with MnSOD, the monooxygenase and cytochrome 
P450 proteins were up-regulated approximately 1.5-fold 
at the end of the exponential phase (Table 1 and addi- 
tional file 4, Fig. S2). These two proteins are closely 
related to the biosynthesis of many secondary metabo- 
lites, including carotenoids [22,23]. Specifically, both cat- 
alyze the addition of a single oxygen atom from 
molecular oxygen to a substrate and the reduction of 
the second oxygen atom into water, a reaction that con- 
sumes two reducing power equivalents. The final donor 
of electrons for the P450 monooxygenases is NADPH 
[44]. Moreover, CrtS (astaxanthin synthase) belongs to 
the cytochrome P450 protein family [45], and CpR has 
recently been identified as an auxiliary enzyme for CrtS 
during astaxanthin synthesis [46]. Two of the proteins 
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identified in this work, cytochrome P450 and monooxy- 
genase, could perform auxiliary reactions during astax- 
anthin biosynthesis; the complete identification and 
further characterization of these proteins is currently 
underway. 

There are clear differences in the induction of astax- 
anthin synthesis between the carotenogenic microorgan- 
isms H. pluvialis and X. dendrorhous. After 24-48 h of 
stress induced by light and high salt, the alga undergoes 
morphological changes and accumulates astaxanthin for 
up to 12 days [43]. In the yeast, under high oxygen con- 
centrations, astaxanthin synthesis is induced on the 
third day of culture, which coincides with the end of the 
exponential phase of growth, and allows the accumula- 
tion of astaxanthin for up to 5 days [22,23]. We found 
similar protein profiles for these microorganisms; how- 
ever, as expected, some of the differentially regulated 
proteins were related to stress response and caroteno- 
genesis. In H. pluvialis, the direct association between 
stress response and carotenogenesis is clear. For X. den- 
drorhous, during aerobic growth with a low level or the 
absence of the antioxidant enzymatic systems, caroteno- 
genesis can be induced. Thus, astaxanthin could per- 
form the antioxidant role of quenching ROS produced 
during cellular metabolism. 

Carotenoid biosynthetic enzymes 

Using our protocol for protein extraction, we deter- 
mined that 9% of all the identified proteins were mem- 
brane associated. We did not identify all of the 
membrane-bound enzymes that perform the late reac- 
tions of carotenogenesis, probably due to technical lim- 
itations. We have identified eight proteins related to 
general or specific steps of astaxanthin biosynthesis. 
Prenyltransferase, geranylgeranyl pyrophosphate 
synthase/polyprenyl synthetase, phytoene desaturase and 
astaxanthin synthase were present similar abundances 
during growth. The other four proteins showed signifi- 
cant fold changes (Table 1 and additional file 4, Fig. 
S2). First, mevalonate kinase (MK), which performs a 
preliminary step in IPP generation, was present at high 
levels in the lag phase, was down regulated in the expo- 
nential phase and was then, up regulated by approxi- 
mately three-fold in the stationary phase (Table 1 
Figure 5 and additional file 4, Fig. S2). This early induc- 
tion is not surprising, as this enzyme performs a preli- 
minary step in common pathways that include 
isoprenoid and ergosterol synthesis. In carotenogenesis, 
it is the second essential enzyme of the mevalonate 
pathway, after 3-hydroxy-3-methyl-glutaryl-CoA reduc- 
tase (HMGR), which catalyzes the phosphorylation of 
mevalonic acid to produce phosphomevalonate. MK 
activity is regulated by intermediates in the pathway, 
such as geranyl pyrophosphate, FPP and GGPP, via 



feedback inhibition [47]. For phosphomevalonate kinase 
we observed the highest abundance at lag phase, while 
diphosphomevalonate decarboxylase reached its highest 
levels during the exponential and stationary phases. 
Because these two proteins perform sequential steps in 
the transformation of mevalonate our results indicate 
that this pathway is tightly regulated to ensure metabo- 
lite availability. 

Another significant carotenoid-synthesis protein is 
phytoene/squalene synthase, which showed higher abun- 
dance at the end of the exponential growth during the 
induction of carotenoid synthesis (Table 1 and addi- 
tional file 4, Fig. S2). This result agrees with our pre- 
viously reported mRNA expression analysis, in which 
the maximal levels of carotenoid-specific genes were 
observed after three days of culture, at the end of the 
exponential growth phase [22,23]. 

In constrast, in H. pluvialis, the mRNA transcript 
levels of carotenoid-related genes reach their maximal 
levels 24-48 h after stress induction, and the synthesis 
and accumulation of astaxanthin occur 6-12 days after 
stress [48]. Another enzyme that performs an initial step 
in carotenogenesis, isopentenyl-diphosphate isomerase 
(IDI), shows maximum expression at 24 h after stress 
induction in H. pluvialis, and is then down-regulated as 
stress persist; a similar behavior has also been observed 
for phytoene desaturase [43,49] (see additional file 3, 
Table S2). Thus, carotenoid-related enzymes in both H. 
pluvialis and X. dendrorhous may have low turnover 
rates; this low rate ensures their long-term activities in 
astaxanthin biosynthesis. 

Conclusions 

In this work, which is the first proteomic characteriza- 
tion of X. dendrorhous, we describe a protocol for the 
enrichment of protein extracts for membrane-bound 
proteins and the efficient extraction of proteins in the 
presence of excess hydrophobic materials such as lipids 
or carotenoids. We have also generated a preliminary 
proteome map, which will be valuable for further studies 
of the organism under different growth conditions. 

We identified two principal types of protein regulation 
associated with astaxanthin biosynthesis. First, we 
observed a tight regulation of carbohydrate and lipid 
proteins, which may ensure the availability of acetyl- 
CoA for the cellular metabolism. Second, we observed 
an up-regulation of redox-specific proteins, such as 
monooxygenase and cytochrome P450, which likely pro- 
vides the redox level necessary for the late reactions of 
astaxanthin synthesis. 

Based on these results, it is possible to assume that 
the production of astaxanthin is an alternative mechan- 
ism for respondings to cellular or environmental stress 
conditions in X. dendrorhous. 
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Although we observed a correlation between mRNA 
levels and protein abundances for phytoene/squalene 
synthase, it will be necessary to perform a membrane 
proteome analysis to study the late enzymes of the 
astaxanthin synthesis. Moreover, detailed transcriptomic, 
proteomic and metabolomic studies are required to gen- 
erate an integrated understanding of the biochemical, 
physiological and biological processes of X. dendrorhous, 
both for basic science research and for metabolic engi- 
neering applications to optimize astaxanthin production. 

Methods 

Preparation of whole-cell protein extracts 

The wild type X. dendrorhous strain ATCC 24230 (UCD 
67-385) was cultured on minimal medium with 2% glu- 
cose as a carbon source [50]. A 10-ml preculture was 
grown to the exponential phase (OD 6.0) at 22°C and 
120 rpm. For the main culture, 250 ml of medium in a 
1-L Erlenmeyer flask were inoculated with 2.5 ml of pre- 
culture and cultivated at 22°C and 120 rpm. For data 
analysis, triplicate cultures in the lag, late exponential 
and stationary growth phases were obtained (Figure 1). 
The cells were harvested by centrifugation at 5,000 x g 
for 10 min at 4°C. After discarding the supernatant, the 
pellet was washed twice with ice-cold water and centri- 
fuged at 5,000 x g for 10 min at 4°C; the washed pellet 
was frozen in liquid nitrogen and stored at -80°C. The 
cell density was determined optically with a spectro- 
photometer at 560 nm and/or gravimetrically by mea- 
suring the cell dry weight. 

Our protein extraction protocol was designed to 
enrich the whole-cell protein extract with membrane- 
bound proteins to allow for the identification of carote- 
nogenic proteins. Yeast cells were lyophilized prior to 
protein extraction. After adding an equal volume around 
500 ul of glass beads (500 urn) to impact-resistant 2-ml 
tubes, the cells were disrupted using a RiboLyzer 
(Hybaid-AGS, Heidelberg, Germany) for 30 s at 4.5 m/s 
and chilled on ice for 1 min to prevent foaming. Five- 
hundred microliters of lysis buffer (100 mM sodium 
bicarbonate, pH 8.8, 0.5% Triton x 100, 1 mM phenyl- 
methylsulfonyl fluoride [PMSF] and protease inhibitors 
[Roche, Mannheim, Germany]) was then added, and the 
samples were incubated for 15 min on ice. Cells were 
disrupted five times for in a RiboLyzer for 30 s at 4.5 
m/s and chilled on ice for 1 min between vortexing 
steps. The cell debris was removed by centrifugation at 
15,000 rpm for 20 min at 4°C, and the supernatant was 
transferred to 1.5-ml tubes. The protein extracts then 
were incubated for 1 h at 4°C with a 10% v/v DNase- 
RNase solution (0.5 M Tris-HCl, pH 7.0, 0.5 M MgCl 2 , 
100 ug/ml RNAse A [Boehringer Mannheim, Germany] 
and 2 ul DNase I [Boehringer Mannheim]). Next, deio- 
nized water was added to produce a final volume of 2.5 



ml, and 200 ul of 0.5 M Tris (pH 6.8) and 20 ul of 1 M 
dithiothreitol (DTT) were added. The samples were 
incubated at room temperature for 30 min. Subse- 
quently, 600 ul of water-saturated phenol was added, 
and the samples were mixed thoroughly and agitated at 
room temperature for 30 minutes. The mixture was cen- 
trifuged at 5,000 rpm at 4°C for 10 min, and the phenol 
phase was transferred into a fresh tube. After the addi- 
tion of 20 ul of 1 M DTT and 30 ul of 8 M ammonium 
acetate, the samples were incubated for 30 min at room 
temperature. The proteins were precipitated by the addi- 
tion of 2 ml of cold (-20°C) methanol and incubation 
over night. The precipitate was centrifuged at 13,000 
rpm at 4°C for 30 min. The supernatant was discarded, 
and the pellet was washed twice with 70% (v/v) cold 
ethanol at -20°C, and incubated for 1 h at 4°C. Finally, 
the pellet was solubilized in 200 ul of buffer (8 M urea, 
2 M thiourea, 2% [w/v] 3[(3-cholamidopropyl)dimethy- 
lammonio]-l-propanesulphonate [CHAPS], 0.01% [w/v] 
bromophenol blue) and stored at -80°C. The protein 
concentration was measured with a Bradford-based pro- 
tein assay (Bio-Rad, Hercules, CA) using bovine serum 
albumin (BSA) as a standard. 

2D electrophoresis 

The resolubilized extract was adjusted to 500 ug in 340 
ul of rehydration buffer, and 1% DTT and 2% immobi- 
lized pH gradient (IPG) buffer at pH 3-10 (IPG buffer, 
Amersham Biosciences, Freiburg, Germany) were added. 
The samples were applied to a 17-cm, non-linear pH 3- 
10 isoelectric focusing (IEF) strip (Immobiline DryStrip, 
Amersham Biosciences) and covered with mineral oil 
(Amersham Biosciences). IEF was carried out on a IPG- 
phor™ system (Amersham Biosciences) using the fol- 
lowing program:10 h at 20°C, 12 h at 30 V, 1 h at 500 
V, 8 h at 1,000 V and 10 h at 8,000 V. 

The strips were equilibrated for 15 min in 10 ml of 
equilibration solution (0.375 M Tris-HCl, pH 8.8, 6 M 
urea, 20% [v/v] glycerol and 2% [w/v] SDS), with 2% (w/ 
v) DTT (reduction step), and for 15 min in 10 ml of the 
equilibration solution with 2% (w/v) iodoacetamide 
(alkylation step). The strip was then applied to a 10% 
SDS-PAGE gel to separate the proteins based on their 
molecular weights (MW). The electrophoresis condi- 
tions were 30 W per gel, applied until the bromophenol 
blue dye front reached the bottom of the gel. 

Protein staining and image analysis 

The gels were fixed in a 10% (v/v) acetic acid and 40% 
(v/v) methanol solution for 2 h, stained for 3 h in a 
Coomassie brilliant blue (CBB) staining solution (2% [w/ 
v] phosphoric acid, 10% [w/v] ammonium sulfate, 5% 
[w/v] CBB G250, 20% [v/v] methanol) and destained 
with 20% (v/v) methanol until the background was clear. 
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The stained gels were scanned and analyzed with 
PDQuest software (version 7.1.1, Bio-Rad). The intensi- 
ties of the spots were normalized to compensate for 
image differences generated by the experimental condi- 
tions (e.g., protein loading or staining) using the total 
density of the valid spots. Spot detection was performed 
using the PDQuest automated spot detection algorithm 
and checked manually. The gel image with the best pro- 
tein pattern and the highest number of spots was cho- 
sen as a reference gel for image analysis, and spots in 
the standard gel were then matched across all gels. To 
compare sets of gels, the MatchSets software tool was 
used to calculate the mean and standard deviation of 
the normalized spot data. For average-fold differences in 
protein abundance, the normalized spot quantity from 
the gel at the lag growth phase was used as a reference; 
the relative abundance levels at later times (i.e., the late 
exponential and stationary phases) were calculated by 
dividing the normalized spot quantity in each gel by the 
abundance data at lag phase. Analyses were validated by 
Student's t-test (p < 0.05). 

MS analyses and database searches 

Coomassie-stained protein spots were excised from the 
2D gels and placed in 96-well plates. The spots were 
destained in 150 ul of 50% acetonitrile (ACN) for 5 min, 
in 150 ul of 50 mM NH4HCO3 and 50% ACN for 30 
min, and then in 150 ul of 10 mM NH4HCO3 for 30 
min while stirring at room temperature. The superna- 
tant was removed, and the plate was dried completely at 
room temperature for 12 h. The proteins were digested 
in-gel with 15 ul of 2.5 mg/ml trypsin (Promega, Madi- 
son, WI) in 10 mM NH4HCO3 at 37°C overnight. 

Samples containing the tryptic peptides were mixed 
1:1 with a solution of 67:33:0.1 water: ACN: trifluoroa- 
cetic acid (TFA) (v/v) saturated with a-cyano-4-hydro- 
xycinnamic acid (CHCA). The mass spectra were 
obtained with an Ultraflex MALDI-TOF-MS (Bruker, 
Bremen, Germany). The spectra data were analyzed in 
detail using FlexAnalysis software (Bruker-Daltonics). 
The peptide mass fingerprints generated by the MALDI- 
TOF MS experiments were interpreted using the Mascot 
search engine run on a local server (Matrix Science, 
London, UK). Each sample was matched to the theoreti- 
cal tryptic digests of proteins from the National Center 
for Biotechnology Information (NCBI) non-redundant 
(nr) database, Swiss-Prot and MSDB. The following 
search parameters were set in the Mascot software: 
taxonomic category, fungi; no MW/pI restrictions; 
enzyme, trypsin; missed cleavages, 1; mass tolerance, 
150 ppm and the modifications of cysteine carbamido- 
methylation and methionine oxidation. The database 
search output contained the number of matched 



proteins ranked according to their Mascot scores, the 
mass error margin and the sequence coverage of the 
matched peptides. A protein was only considered signifi- 
cant if it could be identified at least twice from the same 
position in independent gels, had a Mascot score higher 
than 50 (p < 0.05) and was the same in two of the three 
databases. For the contig identification, the significant 
Mascot score was higher than 40 (p < 0.05). All identi- 
fied proteins were functionally classified according to 
the Kyoto Encyclopedia of Genes and Genomes (KEGG) 
PATHWAY database (http://www.genome.ad.jp/kegg/ 
pathway.html). 

In addition, BLAST (http://blast.ncbi.nlm.nih.gov/ 
Blast.cgi) and CCD conserved domain (http://www.ncbi. 
nlm.nih.gov/Structure/cdd/cdd.shtml) searches were per- 
formed on the predicted or hypothetical proteins that 
had unknown functions to identify structurally and/or 
functionally conserved motifs. 

Carotenoid extraction and HPLC analysis 

Total carotenoids were extracted from the cell pellets 
according to the methods described by An et al. [51]. 
Carotenoids were quantified by absorbance at 465 nm 
with an absorption coefficient of Al% = 2,100. The ana- 
lyses were performed in triplicate, and pigments were 
normalized relative to the dry weight of the yeast. 

Additional material 
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Additional file 1: Fig. SI. 2D gels of soluble proteins from X. 
dendrorhous in the exponential and stationary phases of growth 

Shown are a representative 2D gels for both the exponential and 
stationary growth phases. 

Additional file 2: Table SI. X. dendrorhous proteins identified by 
MALDI-TOF MS. This table lists all MS-identified proteins that were 
separated by 2D electrophoresis. 

Additional file 3: Table S2. Comparative proteomic data from yeast 
and the carotenogenic alga H. pluvialis This table compares the most 
significant results from previous proteomic works on yeast and 
carotenogenic algae. 

Additional file 4: Fig. S2. Differential abundance proteins from X. 
dendrorhous Shown are a representative proteins spots during the 
growth. 
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